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Abstract. Using the semi-classical method of Nikishov-Ritus (NR), the derivation of the transition rate of 
the beams strahlung process is reviewed. This method uses the Furry Picture and the exact solutions of the 
Dirac equation in the external field potential. Por future linear colliders, the nominal machine parameters 
are such that this external field can be considered to be a constant crossed electromagnetic field. The Dirac 
equation solutions can be Pourier transformed such that they are functions of Dirac gamma matrices, Airy 
functions and the usual non-external field solution. The resultant analytic form for the transition rate is the 
same as that obtained by the Quasiclassical Operator (QO) method of Baier-Katkov which is valid in the 
limit of ultra-relativistic electron and vanishingly small radiation angle. The NR calculation however also 
exhibits a pole in the radiation angle for back-radiated photons. The removal of this pole requires a further 
study of IR divergences within the Purry Picture. 



1. Introduction 

There has been a revival in interest over the past decade or so, in fundamental QED processes occur- 
ring in intense external fields. A SLAC experiment involving the interaction of an intense laser with an 
electron beam showed clearly the onset of multiple numbers of laser photons contributing to the process 
Q. Around the same time a similar set of experiments in Europe, on electron laser interaction were 
performed which observed a fermion mass shift due to the action of the electron motion in the laser field 
EO. Another set of experiments has taken place at CERN during the last few years involving fermion 
processes in the intense fields present in the atomic lattices of crystals ||5|. Good agreement has been 
found with theoretical calculations of QED process in crystals Hi. 

The ratio of the strength of the external field to the Schwinger critical field (at which real e+e~ pairs 
can be produced from the vacuum) gives a measure as to whether a quantum or classical treatment of the 
interaction with the external field is necessary. An 'intense field', then, has a field strength sufficiently 
high to require a quantum treatment. Several calculation strategies exist. 

If the field strength is not too high the external field can be considered to contribute one or more 
photons in the Born approximation Q. Higher field strengths require that collective effects be taken 
into account by performing a QED calculation in the Furry Picture (FP). In this picture the external field 
potential is bundled up with the Dirac field to form a bound Dirac field which interacts with the free 
Maxwell field. Exact solutions of the Dirac equation in a (classical) external potential exist and are used 
to calculate the transition probabilities of interest. 



Two calculation strategies are compared in this paper. In the first strategy (also first historically) the 
transition rates of various first order processes in the FP were calculated using the usual S-matrix theory. 
The interpretation of the external field as photons contributing to the process emerges naturally. Exact 
results for transition rates with no kinematic approximations can be obtained. This method is referred to 
as the Nikishov-Ritus (NR) method in this paper ||7l[8l. 

The second method makes use of the fact that in particle colliders the initial state fermions are ultra- 
relativistic and the energy levels of the fermion states in the external field of nominally designed bunches 
are extremely close together. If the fermion final states are also ultra-relativistic then the fermion motion 
can be considered classical in the external field. This Quasi-classical Operator (QO) method starts from 
an intermediate point in perturbation theory before the integration with respect to time has been carried 
out. The transition probability is obtained after allowing the operators of the electron motion to commute 
H. 

The first order processes in the FP for a number of different external fields, including electromagnetic 
fields of various polarisations, magnetic fields and intra-crystalline fields, can be considered well studied. 
Second order processes in the FP require substantial further study, as do the radiative corrections lITOl . In 
this paper the beamsstrahlung process, i.e the bremsstrahlung in the electromagnetic field of a relativistic 
particle bunch, is considered. The NR method for this process will be reviewed and the transition rate 
will be compared to that obtained using the QO method. 

One application of interest for this work is that of spin tracking through the beam-beam interaction 
at colliders. The CAIN program ifTTIl which performs this spin tracking, assumes no radiation angle for 
the outgoing photon. In the interests of precision spin tracking, especially given that collider designs 
include a crossing angle, it is important to know the variation of calculated transition probability with 
radiation angle and no kinematic approximations |[T2l . The full NR result will therefore be examined for 
this behaviour. 



2. QED in intense fields 

The motivation for using the FP arose in the studies of electrons in atomic systems in which bound states 
were natural ones. The interaction of the electrons with the external field can be handled by iteration 
using the Born approximation. However for strong enough fields the Born series does not converge and 
the interaction with the external field has to be taken into account to all orders. The FP achieves this 
by repackaging the Lagrangian density of the system with the interaction term specifying the coupling 
between free Maxwell and the combined Dirac-external field system. The usual perturbation method 
within S-matrix theory is then employed using the solution for the bound Dirac field equation. 
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The solution i/jy{x,p) to the bound Dirac field equation was worked out long ago |[T3ll . A nice 
exposition of its derivation is found in section 40 of |[T4l . This Volkov solution is a product of a Volkov 
Ep{x) function of the external field 4-potential A^, and the usual free field solution. 
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The orthogonality of the Volkov Ep{x) is an important result |[l5l[T6l which allows the bound fermion 
propagator x') to be written as the usual propagator sandwiched between Volkov Ep{x) functions, 
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Using the above results the required interaction between bound fermions and photons can be studied 
using the usual S-matrix theory and Wick's theorem. Since the Volkov solution is written in terms of the 
bispinor and gamma matrices the usual Dirac gamma algebra and spin and polarization sums are used to 
obtain the transition probability of the desired QED process. This was the route chosen by Nikishov and 
Ritus (and others) in the 1960s and which is summarised in the following two sections. 

With the solutions of the bound Dirac equation and an expression for the bound fermion propagator, 
Volkov Ep{x) functions can be grouped together to form only one new element in the S-matrix - a 
modified vertex 

= y d^x Epf {x)-f^Ep^ (x) exp [z(p^ -Pf - kf) • x)] (4) 

In terms of the Feynman diagrams, the bound Dirac equation solutions are represented as double 
straight lines. 

3. Bremsstrahlung in a constant crossed field 

The beamsstrahlung process (figure [T]) is simply the bremsstrahlung process that takes place when a 
charge radiates a photon in the field of an oncoming charge bunch in collider collisions. For proposed 
future linear colliders the radiation coherence length - classically, the ratio of the radius of curvature 
of the radiating charge to the relativistic gamma value of the radiating charge - is shorter than both the 
typical bunch length and the characteristic length specified by the beam disruption parameter. Therefore 
the external field can be considered to be a constant crossed field ifTTIl specified by 
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The Lorentz condition and a choice of reference frame such that the zeroth (the time) component 
of the external field 4-potential is zero, requires that the 3-potential ai and the 3-momentum k of the 
external field be orthogonal. Using these vectors as a basis for a coordinate system with di forming the 
X-axis and k the z-axis, the radiation angles {Of^cpf) of the radiated photon kf can be specified (figure 
[2]). If the radiating electron is directed exactly opposite to the wave- vector of the oncoming bunch field 
(as is usually the case in colliders except for a small crossing angle) then the angle Of = n corresponds 
to forward radiation and = to back radiation 

An explicit form for the Volkov Ep(x) functions can be obtained by substitution of the 4-potential of 
the constant crossed field. In order to simplify the dependence on the space-time variable x^, a Fourier 
transform is made to an infinite integration over a variable r 



Figure 1. Beamsstrahlung Feynman 
diagram. 



Figure 2. Coordinate system and radia- 
tion angles. 
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The modified vertex can Hkewise be transformed and the integration over space-time carried out to 
produce the usual delta function whose argument expresses the momentum conservation 



There is a contribution rk^ from the external field, which is integrated over all values allowed by con- 
servation of momentum. It is open to interpretation whether this is indeed one photon of 4-momentum 
rkfj^ or r photons of 4-momenta A;^ each. This is unlike the case for processes in laser fields in which 
the contribution from the external field is a summation (over, say, n) rather than an integration and 
the interpretation clearly is in terms of n external field photons. Though external field photons were not 
included at the outset, they emerge through the quantum interaction of the electron with the external field. 

2 2 

For a constant crossed field, the concept of a quasi-momentum for the fermions g'/i = + ^/i 

2 2 

is ill defined, since the intensity of the crossed field oo. Consequently the conservation of 

momentum of the process is written in terms of the usual external field- free momentum p^. 

4. Beamsstrahlung transition probability 

The matrix element of the beamsstrahlung process is the modified vertex sandwiched between Dirac 
bispinors. In deriving the differential transition probability, dw, it proved useful to use relativistic 
normalization for the field operators and to follow the derivation of the general form in ifTSl . It should 
be noted that the square of the matrix element will introduce another integration (over variable r') and a 
second delta function. The particle energy-momenta are specified by (ci^pi), (cf^pf) and {uf.kf). 
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There is nothing conceptually difficult in reducing the transition probability to its simplest form, 
however the procedure is technical and time-consuming. The main steps are simply listed and the end 
result quoted for final analysis: 

• Perform the integration over initial 3-momenta using 3 components of one of the delta functions 

• Write the Fourier transformed functions F„ as Airy functions by using the relation 
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• Perform the spin sums and trace calculation using the usual Dirac algebra 

• Introduce an additional delta function using the relation 



2e/ 

• Using two delta functions do the integrations over final fermion 3-momenta and the contribution 
from the external field r, to get the condition r ^ {pi • kf)/{k • pf) 

• Integrate over the x-component of the final photon momentum to get the condition that r' ^ r 



and transform the transition probability it; to a transition rate W 



• Transform the y-component of the final photon 3 -momentum fc^ and integrate to reduce products of 
Airy functions to a single Airy function |[T9l using 

/ ^Ar (t + a) = - / Ai(t) dt 
Jo Vt 2J22/3a 

• Replace Airy functions with Bessel functions using relations like, Ai(t) = -^^1/3 (1*^^^) 

• The integration over the z-component of the final photon 3-momenta is transformed to a function 
of scalar products of 4-momenta u = ^^'^^j 



The final result for the transition rate is 

^2 1 



dW 
du 



am 



7r\/3ei (l + n)2 
2u 



l + u 



^2/3 (x) 



(9) 



5. Comparison of the Nikishov-Ritus and Baier-Katkov methods 

The same beamsstrahlung transition rate appears in the Uterature for both NR and QO methods, though 
the exact form of the integration variable varies from derivation to derivation and NR prefers to write in 
terms of Airy functions. Nevertheless, a direct comparison can be made and the difference in the final 
result for the two methods is encapsulated in the integration variable 

x(QO) ^ ^ 
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Physically, the parameter x is inversely proportional to the magnetic field strength of the oncoming 
charge bunch and proportional to the ratio of radiated photon energy to final electron energy. The x(NR) 
also retains an angle dependence in its scalar products. Naively, if the electron is ultra-relativistic and 
radiates almost entirely in the forward direction 0j = tt (as one might expect) then x(NR) x(QO) 
and the two results seem consistent in the limit assumed by the QO method. 

However, the beamsstrahlung transition rate contains a pole at the lower end of the x integration (at 
X = 0) since the K Bessel function of fractional order approaches infinity as its argument approaches 
zero i20|. Comparison between the QO and NR methods near the pole requires careful consideration. 
Both x(NR) and x(QO) t^i^d towards zero as the energy of the radiated photon becomes small. This is 
the IR divergence. Applying the usual argument for the removal of the divergence, a cut is made at a 
minimum energy which is the limit of the energy resolution of the detector. The residual divergent pole 
is then removed theoretically by inclusion of radiative corrections, though it should be noted that it does 
not appear to have been done explicitly in the literature for this process in the FR 

Writing the minimum photon energy cut instead as a fictions small photon mass the parameter 
x(NR) can be written in terms of the radiation angle 
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and its dependence plotted (figure [3(a)] ). The closer x(NR) reaches zero, the more closely the transi- 
tion rate pole is reached. Clearly this is the case for 9f = 0, i.e. back radiation and is the case even for 
ultra-relativistic electrons. This is not a trivial observation since even with an energy cut, the transition 
rate reaches a maximum around the pole and any small variation at that point will have a relatively large 
effect on the transition rate value. 



A maximum transition rate for back radiation is obviously in contradiction to the results of classical 
electrodynamics which predict forward radiation for ultra relativistic electrons. Additionally the pole in 
the transition rate for back radiation renders the equivalence of the QO and NR results within the ultra 
relativistic limit invalid. 



In principle, the removal of the IR divergence from the beamsstrahlung process within the FP should 
be straightforward. It has been argued that since IR divergences cancel in the normal interaction picture, 
then it must also be the case for the bound interaction picture ll22l . The removal may require a calculation 
of the vertex function in the FR In any case a dedicated further study is necessary. 
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Figure 3. Differential transition rate 



dW 
dx(NR) 



and parameter x(NR) near the transition rate pole. 



6. Conclusion 

The analytic form of the transition rate of the beamsstrahlung process has been compared for two calcula- 
tion methods. The QO method assumes a classical interaction between the electron and the external field, 
but a quantum interaction to produce the radiated photon. The spectrum of the radiation (the derivative 
of the transition rate with respect to the radiated photon energy) shows the usual IR divergence for low 
energy photons. The QO calculation is only expected to be valid for an ultra-relativistic electron before 
and after the radiation. 



The NR method is a full calculation using the bound interaction picture within the usual S-matrix the- 
ory with no kinematic approximations. In certain limits the analytic form of the NR radiation intensity 
appears to reduce to that of the classical Schott formula ll2ll and to the beamsstrahlung transition rate 
obtained by the QO method. However a full consideration of the NR calculation reveals a transition rate 
maximum for back-radiated photons. It seems clear that an explicit inclusion of radiative corrections 
within the FP is required and such a work is intended for a forthcoming paper. 
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